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ABSTRACT 

A new method of matrix template matching is presented in the context of pulsar timing analysis. 
Pulse arrival times are typically measured using only the observed total intensity light curve. The new 
technique exploits the additional timing information available in the polarization of the pulsar signal by 
modeling the transformation between two polarized light curves in the Fourier domain. For a number 
of millisecond pulsars, arrival time estimates derived from polarimetric data are predicted to exhibit 
greater precision and accuracy than those derived from the total intensity alone. Furthermore, the 
transformation matrix produced during template matching may be used to calibrate observations of 
other point sources. 



Subject headings: methods: data analysis 
polarimetric 

1. INTRODUCTION 



polarization — pulsars: general — techniques: 



High-precision pulsar timing is a well-established tech- 
nique of modern astrophysics; it has yielded the strongest 
constraints on theories of gravitation in the strong-field 
regime (Stairs 2004), and is anticipated to provide a direct 
detection of the stochastic gravitational wave background 
due to supermassive black hole binaries (Jenet et al. 2005). 
Fundamental to every pulsar timing experiment is a mea- 
surement known as the pulse time-of-arrival (TOA), the 
epoch at which a fiducial phase of the pulsar's periodic sig- 
nal is received at the observatory. The confidence limits of 
the various physical parameters that are derived from the 
TOA data depend upon the precision and accuracy with 
which arrival times can be estimated. 

In addition to typical constraints such as the system 
temperature, instrumental bandwidth, and integration 
length, TOA precision also depends upon the physical 
properties of the observed pulsar, including its total flux 
density, pulse period, and the shape of its phase-resolved 
light curve, or pulse profile. When fully resolved, nar- 
row features in the pulse profile provide strong constraints 
during template matching, thereby yielding greater ar- 
rival time precision (see § 3.2). The polarized component 
of the pulsar signal often displays much sharper features 
than observed in the total intensity, especially when the 
pulsar exhibits transitions between orthogonally polarized 
modes. This important property may be exploited to sig- 
nificantly improve timing precision by incorporating po- 
larization data. 

Pulsar polarization also impacts on timing accuracy 
(Cordes et al. 2004). The total intensity profile can be 
significantly distorted by instrumental artifacts, a problem 
most readily observed as a systematic variation of arrival 
time residuals with parallactic angle (van Straten 2003). 
To address this issue, Britton (2000) proposed timing the 
polarimetric invariant profile, which greatly improved the 
timing accuracy of PSRJ0437-4715 (Britton et al. 2000). 
However, there are disadvantages of using the invariant 
profile; its signal-to-noise ratio (S/N) is at most 1/a/2 



times that of the total intensity (see § 4) and inversely 
proportional to the degree of polarization; also, its com- 
putation suffers from imprecision in the estimation of the 
off-pulse signal. Therefore, use of the invariant profile can 
be detrimental when the pulsar is highly polarized or has 
low flux density. 

In contrast, the matrix template matching technique 
presented in this paper can be used to improve both the 
precision and accuracy of arrival time estimates with- 
out any sensitivity to mean off-pulse polarization. The 
method simultaneously yields the polarimetric transfor- 
mation between the template and observed profiles, which 
may be utilized to completely calibrate the instrumental 
response in other observations. Following a review of the 
required mathematics in § 2, the matrix template matching 
technique is formulated and quantitatively compared with 
scalar methods in § 3. In § 4, the analysis is demonstrated 
using a sample of millisecond pulsars. The application of 
matrix template matching to instrumental calibration is 
outlined in § 5, and the main conclusions of this work are 
summarized in § 6. 

2. REVIEW OF THE JONES CALCULUS 

The formulation and analysis of the matrix template 
matching method utilizes the terminology and notation 
reviewed in this section. The polarization of electromag- 
netic radiation is described by the second-order statistics 
of the transverse electric field vector, e, as represented us- 
ing the complex 2x2 coherency matrix, p = (e (g) e^) 
(Born & Wolf 1980). Here, the angular brackets denote 
an ensemble average, (g is the matrix direct product, and 

is the Hermitian transpose of e. The coherency matrix 
is commonly expressed as a linear combination of Hermi- 
tian basis matrices, p = {Sq ctq + S • t)/2, where (Tq is the 
2x2 identity matrix, cr — (cri, cr2, cr^) are the Pauli spin 
matrices, Sq is the total intensity, and S ~ {Si, S2, S3) is 
the polarization vector (Britton 2000) . The Pauli matrices 
are traceless and satisfy cr^ = ctq; therefore, Sk = tr((Tfep), 
where tr is the matrix trace operator (Hamaker 2000). 

In the analysis of the reception of polarized radiation, 
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the following conventions are used. The response of a sin- 
gle receptor is defined by the Jones vector, r, such that 
the voltage induced in the receptor by the incident elec- 
tric field is given by the scalar product, v = r^e. A dual- 
receptor fced is represented by the Hcrmitian transpose of 
a Jones matrix with columns equal to the Jones vector of 
each receptor. 



{ro ri)^ 



' 00 

'■1*0 



(1) 



The receptors in an ideal feed respond to orthogonal senses 
of polarization (ie. the scalar product, r^ri = 0) and have 

identical gains (ie. rlro = r\ri). 

A more meaningful geometric interpretation of Jones 
vectors is provided by the corresponding Stokes parame- 
ters, Sk = tr[(Tfer r^]. The state of polarization to which 
a receptor maximally responds is completely described by 
the three components of its associated Stokes polarization 
vector, S. Therefore, it is common to define a receptor us- 
ing the spherical coordinates of S (Chandrasekhar 1960); 
in the linear basis, these include the gain, g = \r\ = \S\^, 
the orientation. 



and the ellipticity. 



such that 



n 1 , _i 5'2 
^=-tan -, 



1 . -1 ^3 



COS 6 cos e -|- i sin ^ sin e 

sin 9 cos e — i cos 6 sin e 



(2) 



(3) 



(4) 



The impact of non-ideal feed receptors on pulsar timing 
is analyzed by exploiting a powerful classification of Jones 
matrices motivated by the polar decomposition. Any non- 
singular matrix can be decomposed into the product of a 
unitary matrix and a positive-definite Hermitian matrix. 
Using the axis-angle parameterization (Britton 2000), the 
polar decomposition of a Jones matrix (Hamaker 2000) is 
expressed as 

J = JB^(/3)Rft(<^), (5) 

where J = (det J)5, B^(/3) is positive-definite Hermitian, 
and Rft((/)) is unitary; both B.^(/3) and Rft(0) are uni- 
modular. The unit 3-vectors, rh and n, correspond to axes 
of symmetry in the three-dimensional space of the Stokes 
polarization vector. Under the congruence transformation 
of the coherency matrix, the Hermitian matrices, 

Brh(/3) = exp(/3 m • cr) = aQ cosh/3 + rh • a sinh (3, (6) 

effect a Lorentz boost of the Stokes 4- vector along the rh 
axis by an impact parameter 2/3, such that the resulting 
total intensity, Sq = So cosh 2(3 + S • rh sinh 2/3. Similarly, 
the unitary matrices, 

R-Ti(0) = exp{i(j)h • a) = ao cos(t) + in • <Tsin</), (7) 

rotate the Stokes polarization vector about the n axis 
by an angle 2(f>, leaving the total intensity unchanged. 
This parameterization enables the important distinction 
between polarimctric transformations that mix the total 
and polarized intensities (boosts) and those that effect a 
change of basis (rotations). 



3. MATRIX TEMPLATE MATCHING 

The primary purpose of this paper is the formal de- 
scription of the matrix template matching technique and 
the quantitative comparison of its effectiveness with that 
of conventional scalar methods. The performance of a 
method of TOA measurement may be evaluated in terms 
of the precision and accuracy of the arrival time estimates 
that it produces. The analysis of TOA precision requires 
careful attention to the propagation of experimental er- 
ror, as described in § 3.2. TOA accuracy depends upon 
the susceptibility of the technique to sources of system- 
atic error. To quantitatively compare different methods, a 
simulation that spans the full range of potential artifacts 
is devised in § 3.3. 



3.1. Description of Technique 

A pulsar's mean pulse profile is measured by integrat- 
ing the observed flux density as a function of pulse phase. 
By averaging many pulse profiles, one with high S/N may 
be formed and used as a template against which the indi- 
vidual observations are matched. The best-fit phase shift 
derived by the template matching procedure is then used 
to compute the pulse TOA. 

Taylor (1992) presents a method for modeling the phase 
shift between the template and observed total intensity 
profiles in the Fourier domain. In the current treatment, 
the scalar equation that relates two total intensity pro- 
files is replaced by an analogous matrix equation, which 
is expressed using the Jones calculus. Let the coherency 
matrices, p'(0„), represent the observed polarization as a 
function of discrete pulse phase, where < n < N 
and N is the number of intervals into which the pulse pe- 
riod is evenly divided. Each observed polarization profile 
is related to the template, Poi^n), by the matrix equation. 



p'(0„) = Jpo( 



^)jt 



Pdc + PnI 



0, 



(8) 



where J is the polarimctric transformation, ip is the phase 
shift, pj^Q is the DC offset between the two profiles, and 
represents the system noise. The discrete Fourier trans- 
form (DFT) of equation (8) is 

p'(t'm) = JPoi'^m)^^ exp{-i2TrUmV>) + P^ivm), (9) 

where is the discrete frequency. Given the observed 
Stokes parameters, S'i.{(j)n), and their DFTs, S'^(t'm)5 the 
best-fit model parameters will minimize the objective 
merit function. 



N/2 3 
m=l k=0 



(10) 



where <;k is equal to the rms of the noise in each DFT and 
tr is the matrix trace operator. As in van Straten (2004), 
the partial derivatives of equation (10) are computed with 
respect to both ip and the seven non-degenerate parame- 
ters that determine J. The Levenberg-Marquardt method 
is then applied to find the parameters that minimize 
(Press et al. 1992). 
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3.2. Timing Precision 

When compared with the scalar technique, matrix tem- 
plate matching quadruples the number of observational 
constraints while introducing only six degrees of freedom. 
Therefore, arrival time estimates derived from the polar- 
ization profile might be expected to have greater precision 
than those derived from the total intensity profile alone. 
However, the effectiveness of matrix template matching de- 
pends upon both the degree of polarization and the vari- 
ability of the polarization vector as a function of pulse 
phase. These properties determine the extent to which 
the phase shift, (p, is correlated with the free parameters 
that determine the Jones matrix, J. 

To estimate the timing precision attainable by matrix 
template matching, consider the solution of equation (8) 
in the special case that J is known. After calibration, the 
minimization of equation (10) by variation of (p requires 
finding the appropriate root of 



0)1 
dip 



47r ^ ^ ''j^ Vm sin((/)fc,, 

m=l fe=0 



^1 



2nipi/m), (11) 



where Sk,m = 'S'^*(^'m)tr[crfe Po{i^m)] is the cross-spectral 
power of the template and observation and (l)k,m is the 
argument of Sk,m- To first order, 



Q 2 N/2 3 



^k,m I 



i(t>k,mT^m - = 0, (12) 



which is equivalent to the solution of a line with slope 

2Trip passing through the origin and the points, {i>m, 4'k,m)- 
This linear approximation is readily solved for tp and used 
to determine the conditional variance, 



var((p|J) = 



N/2 



m=l 



k=0 



(13) 



Equation 13 demonstrates that fluctuation power con- 
tributes quadratically as a function of frequency to the 
reduction of the phase shift variance. In other words, 
sharper features in the pulse profile, which generate more 
power at higher harmonics, yield greater arrival time preci- 
sion. This general results holds for both matrix and scalar 
template matching (in the scalar case, the sum over index 
k stops at 0). 

Equation 13 also provides an upper limit on the arrival 
time precision that may be obtained in the ideal case of a 
perfectly calibrated instrument. However, the instrumen- 
tal response is generally unknown, and the Jones matrix 
must also be varied in order to minimize equation (10). 
The covariances between the phase shift, and the free 
parameters that describe J will increase the uncertainty in 
ip and therefore decrease arrival time precision. 

Formally, the variance of is determined by the covari- 
ance matrix of the free parameters, C = cx~^, where a is 
the curvature matrix. Assuming that the noise power in 
the DFT of each Stokes parameter is equal, 

N/2 



1 d'x' 



2 dr]rdr]s 



tr 



dr]r dr]s 



(14) 



If the free parameters, 77, are partitioned into (f = rjo and 
the seven Jones matrix parameters, rjj, then C may be 



conformably partitioned into (Stuart, Ord, & Arnold 1999) 
C = ( ^^^^ ) , (15) 

and yaT{ip) = c^^. Furthermore, if the free parameters are 
multinormally distributed, then the multiple correlation 
between (p and 77 j. 



(16) 



describes the relationship between the variance and con- 
ditional variance of (p, 

var(^|J)=var(^)(l-i?2j). (17) 

The multiple correlation coeSicient, < R^pj < 1, provides 
a useful measure of the decrease in arrival time precision 
that results from an unknown instrumental response. It is 
important to note that R^j can be computed using only 
the template polarization profile; that is, though a de- 
pends on and c, all factors of <; are canceled in equa- 
tion (16). Therefore, R^j is unique to each pulsar and is 
independent of the S/N of the observations. 

However, the precision and accuracy with which R^pj 
can be estimated does depend on the S/N of the template. 
In fact, measurement noise in the template polarization 
profile artificially decreases both the multiple correlation 
and the conditional variance of ip. To minimize the im- 
pact of noise on the computation of var((p|J), R^j, and 
var((p), the summations over the index m in equations (13) 
and (14) are performed up to a maximum harmonic, M, 
the highest frequency at which the fluctuation power spec- 
tra exhibit three consecutive harmonics with power greater 
than three times the mean noise power. 

To completely eliminate the dependence on var(<^|J) 
and var((^) are normalized by the corresponding variance 
in the phase shift, (p, yielded by the conventional method of 
scalar template matching the total intensity profile. That 
is, equation (13) is used to predict the relative conditional 
timing error in the case of a known instrumental response. 



(18) 



where G is the absolute gain, 

M 



var((p J) 


2 


[ So 1 


var{(p\G)_ 




So -|- Si_3_ 



So = ^ i^llSoM' 



and 



Similarly, 



^1-3 



M 



Cm = 



(19) 



(20) 



(21) 



yields the relative timing error in the general case of an 
unknown response. Given a template polarization pro- 
file, equations (13) - (21) can be used to predict arrival 
time uncertainty in the case of either known or unknown 
polarimetric response. These equations provide the basis 
for comparing the precision of scalar and matrix template 
matching methods for each pulsar in § 4. 

Note that (T^|j < 1; that is, for a known response, the 
precision yielded by matrix template matching will always 
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be as good as or bettor than that of the conventional scalar 
method. In the general case of unknown response, matrix 
template matching will yield more precise TOA estimates 
only when det ex ^ and 

^'j+<^'|j<l- (22) 

The interpretation of the first restriction is simplified by 
considering the polar decomposition of the unknown re- 
sponse, J. To constrain the rotation component of this 
matrix, the pulse profile must include at least two non- 
collinear polarization vectors. Therefore, at becomes sin- 
gular when either the profile is unpolarized or the polar- 
ization vectors. S{4>n), lie along a single line. In these 
two special cases, it is not possible to invert a and matrix 
template matching fails. 

Equation (22) places an upper limit on the value of 
Rip J . As the degree of polarization of the pulse profile ap- 
proaches zero, (T^|j approaches unity, and the maximum 
value of Ripj permitted by equation (22) approaches zero. 
Owing to the unknown boost component of the response, 
the multiple correlation approaches unity when there is a 
high degree of symmetry in S'o (</>«) combined with anti- 
symmetry in S{(j)n). Similarly, axial symmetry in S{(j)n) 
will increase the multiple correlation between the rotation 
component and <p. In these cases, the uncertainty of ar- 
rival time estimates derived from the polarization profile 
will be larger than that of those derived from the total 
intensity alone. 

In summary, matrix template matching does not per- 
form as well as the conventional method when the degree 
of polarization is low and the multiple correlation (due 
to the symmetry properties of the polarization profile) 
is high. However, in general, the radiation from pulsars 
is highly polarized and, especially in the millisecond pul- 
sar population, most polarization profiles exhibit complex 
structure. Therefore, it is expected that, in all but a few 
exceptional cases, matrix template matching will produce 
arrival time estimates with greater precision than those 
derived by conventional scalar methods. 

3.3. Timing Accuracy 

Unmodeled polarimetric distortion of the total inten- 
sity profile can shift arrival time estimates derived by con- 
ventional methods. As equation (8) incorporates an un- 
known polarimetric transformation, it has the capacity to 
model instrumental polarization instabilities and isolate 
them from pulse arrival time variations, thereby eliminat- 
ing a potentially significant source of systematic timing er- 
ror. The improvement in timing accuracy may be analyzed 
by designing a simulation in which arrival time estimates 
are derived from polarization profiles that have been sub- 
jected to transformations that distort the total intensity 
profile. 

As reviewed in § 2, only boost transformations affect 

the total intcmsity. Physically, boosts arise from the dif- 
ferential amplification and non-orthogonality of the feed 
receptors (Britton 2000, and references therein). A brief 
review of these phenomena provides the motivation for the 
analysis technique presented in this section. 

For a pair of orthogonal receptors with different gains, 
go and gi, define the orthonormal receptors, ro = fo/go 
^ There is an error following equation (14) in Britton (2000), where it 



and fi — ri/gi, and substitute into equation (1) to yield 
J = (ro ri)t = G I ) (rori)t, (23) 

where G = (,go.9i)^ is the absolute gain and F = (.go/.9i)^ 
parameterizes the differential gain matrix. Equation (23) 
is a polar decomposition and, by substituting F = exp(/J), 
the differential gain matrix may be expressed in the form 
of equation (6), 

B^(/3) = ( )^rTo cosh/3 + sinh/3, (24) 



where rh = (1, 0, 0) and^ 

/3 = 



lln^. 

2 gi 



(25) 



In the case of linearly polarized receptors, the rh axis lies 
in the Stokes Q-U plane; for circularly polarized receptors, 

rh corresponds to ± Stokes V . To first order, [3 = 7/2, 
where 7 = ffo/Si ~ 1 is the differential gain ratio. 

For a pair of non-orthogonal receptors, first consider the 
spherical coordinate system introduced in § 2. The ori- 
entations and ellipticities of orthogonal receptors satisfy 
do — di = ±7r/2 and eo = — ei. If 5g and (5^ parameterize 
the departure from orthogonality in each of these angles, 
such that ^0 = ^0 + and eg = eo -|- ^e, then 

TqVi = gogi{smSg cos{Se + 2eo) + icosSgtiinS^). (26) 
However, this description is specific to linearly polarized 
receptors. In the circular basis, eo =^ 7''/4, and the non- 
orthogonality is completely described by 5^; Sg becomes 
degenerate with the differential phase of the receptors. As 
such a degeneracy exists at the poles of any spherical co- 
ordinate system, a basis-independent parameterization of 
non-orthogonality is sought. 

To this end, it proves useful to consider the geometric 
relationship between the Stokes polarization vector of each 
receptor. In particular, the scalar product, 

So-Si=2|rJri|2-|ro|2|ri|2, (27) 
shows that orthogonally polarized receptors have anti- 
parallel Stokes polarization vectors {Sq-Si = — |So||<S'i|). 
Furthermore, where 8 is the angle between Sq and S'l, 
the angle S = {n — 6)/2 parameterizes the magnitude of 
the receptor non-orthogonality, such that 



sin (5 



\rlri\ 



(28) 



|T-o||ri| 

It is much simpler to relate 6 to the boost transformation 
that results from non-orthogonal receptors. To determine 
the boost component of an arbitrary matrix, J, the po- 
lar decomposition (eq. [5]) is multiplied by its Hermitian 
transpose to yield 

JJ^i = |detJ|B^(/3) = I det J I Brf, (2/3). (29) 
For a pair of receptors with gain, g, substitution of equa- 
tion (1) into equation (29) yields 

"n, (30) 



Jjt - ^ 



9 



where W = rlvi. Substitute W ■ 



2 

: g^e 



' tanh 2/3, so that 



I det J I = det(JJt)5 = (5^ _ \w\'^)i = g2sech2l3, and 
jjt 

B^(2/3) = J- — =-r = CTo cosh 2/3 rh • crsinh 2/3, (31) 

should read p = ln(pa/56)/2. 
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where rh = (O,cos$,sin<i>) and 

/3=^tanh-i^^. (32) 
2 

To first order, /3 ~ 5/2 (see eq. [28]), which is consistent 
with the approximation in equation (19) of Britton (2000). 

Due to the combined effects of differential gain and re- 
ceptor non-orthogonality, the boost axis, m, can have an 
arbitrary orientation. For example, to first order in the 
polar coordinate system best-suited to the linear basis, 
rh Gc (7, (5e, iJe). Furthermore, the instrumental boost can 
vary as a function of both time and frequency for a va- 
riety of reasons. For example, the parallactic rotation of 
the receiver feed during transit of the source changes the 
orientation of rh with respect to the equatorial coordinate 
system. Also, to keep the signal power within operating 
limits, some instriiments employ active attenuators that 
introduce differential gain fluctuations on short timescales. 
Furthermore, the mismatched responses of the filters used 
in downconversion typically lead to variation of 7 as a 
function of frequency (see § 5). 

The impact of these variations on conventional timing 
accuracy may be estimated through a simulation in which 
copies of the template polarization profile are subjected to 
a boost transformation before the phase shift, (f, between 
the template and the distorted copy is measured (using 
only the total intensity). By varying the orientation of the 
boost axis, m, the maximum and minimum phase shift off- 
sets for a given level of distortion, /3, are found and used 
to define the systematic timing error, 

Ar(/3) = P(max{^(/3, rh)} - min{^(/3, rh)}), (33) 
where P is the pulsar spin period. Empirical observations 
made during testing show that, to first order, At varies 
linearly with /3. 

4. APPLICATION TO SELECTED PULSARS 

To demonstrate the results of the previous sections, 
the analysis is applied to a sample of millisecond pul- 
sars among the best for high-precision timing experiments. 
Polarization profiles with high S/N were obtained for 
PSR J0437-4715 (van Straten 2005), PSR J1022-I-1001, 
PSR J1713-h0747, PSR J1909-3744 (Ord et al. 2004), 
PSR B1855-I-09 (I. Stairs 2005, private commimication) 
and PSR B1937-I-21 (Stairs, Thorsett, & Camilo 1999). 
For each pulsar, the fluctuation power spectra of the total 
intensity, \So{u,„)\'^ , and polarization, |S'(i'„i)p, are plot- 
ted in Figure 1. The maximum harmonic used to com- 
pute the arrival time uncertainty for each pulsar (as de- 
scribed in § 3.2) is indicated by a vertical line in each panel. 
Over some frequency intervals, the polarization fluctuation 
power of PSR J0437-4715 and PSR J1022-^-1001 exceeds 
that of the total intensity. Therefore, it is expected that 
out of the selected sample these two pulsars will benefit 
most from the application of matrix template matching. 

The values of a^,, R^j, and a^^j (see § 3.2) for each of the 
selected pulsars are listed in Table 1. In all cases, arrival 
time precision is predicted to improve through matrix tem- 
plate matching, with timing errors decreasing between 4% 
and 33%. As expected, the pulsars with the greatest im- 
provements (PSRJ0437-4715 and PSR J1022+1001) are 
also those with the greatest amount of polarization fluctu- 
ation power (relative to that of the total intensity) at high 
frequencies. 



For further comparison, Table 1 also includes a^p, the 
relative uncertainty of the arrival time estimate derived 
by scalar template matching the invariant profile, 

SinAM = {[SoiM? - Wn)]^) ^ . (34) 

Again, is normalized with respect to the uncertainty of 
the total intensity TOA. As predicted by a simple consid- 
eration of the noise power in Sinv{4'n), the uncertainties 
of the invariant TOAs are at least \/2 times greater than 
those of the total intensity TOAs. For PSR .11022+1001, 
the precision yielded by matrix template matching is al- 
most three times better than that of the invariant TOA. 

The predicted values of in Table 1 represent the 
minimum increase in experimental sensitivity yielded by 
matrix template matching. They do not include the po- 
tentially significant improvements in timing accuracy that 
may be gained by using this technique. For example, the 
difference between predicted value for PSR J0437— 4715, 
= 0.8318(1), and the measured value of ~ 0.74 reported 
in van Straten (2005) can be explained by the presence 
of systematic timing error in the conventionally-derived 
TOAs of that experiment. 

Timing accuracy is addressed in Table 2, which lists the 
best published rms timing residual, Cr, for each pulsar and 
the magnitude of instrumental artifacts that will produce 
significant systematic timing error in each data set, as de- 
fined by At ~ ar (see eq. [33]). Note that systematic 
timing errors exist only in the TOAs derived from the to- 
tal intensity profile. The simulation confirms that, within 
the experimental uncertainty, there is no distortion of ar- 
rival times derived using matrix template matching, the 
free transformation J completely absorbs the simulated 
instrumental boost. 

In Table 2, the magnitude of the instrumental error is 
specified using the boost impact parameter, Ps, the equiv- 
alent differential gain ratio, 7^ = 2(3s, and the equiva- 
lent receptor non-orthogonality, Ss = 2f3s- Except for 
PSR J1713-F0747 and PSR B1855+09, a conservative 
value of 7s ^ 3% or Sg ^ 2° results in significant sys- 
tematic distortion of arrival times derived by conventional 
methods. Not surprisingly, the two pulsars predicted to be 
most susceptible to instrumental error (PSR J0437— 4715 
and PSR J1022+1001) are also those for which the sys- 
tematic distortion of arrival time estimates due to poor 
calibration has already been reported (van Straten 2003; 
Hotan, Bailes, & Ord 2004). 



5. INSTRUMENTAL CALIBRATION 

In addition to the phase shift, (p, the matrix template 
matching method simultaneously yields the polarimetric 
response, J, required to transform the observation into 
the template. Assuming that the template has been well- 
calibrated, this unique transformation may be applied 
to completely calibrate the instrumental response in ob- 
servations of other point sources. This is demonstrated 
by matching an uncalibrated, five minute integration of 
PSR J0437— 4715 to its standard template, producing the 
instrumental parameters shown with their formal standard 
deviations in Figure 2. As in van Straten (2004), the in- 
strumental response is parameterized by the orientations. 



6 



W. van Straten 




Fig. 1. — Fluctuation power spectra of total intensity (solid) and polarization (dotted). The vertical line that extends half the height of 
each plot indicates the maximum harmonic used to estimate the arrival time precision. The same number of harmonics are plotted for each 
pulsar, except PSR J0437— 4715 and PSR J 1909-3744, for which twice as many harmonics are required. 
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Table 1 

Relative Arrival Time Uncertainties 



Pulsar 










J0437-4715 


0.8285(1) 


0.08993(8) 


0.8318(1) 


1.4314(3) 


J1022+1001 


0.615(4) 


0.394(1) 


0.669(4) 


1.80(2) 


J1713+0747 


0.877(2) 


0.047(3) 


0.877(2) 


1.543(5) 


B1855+09 


0.9263(9) 


0.1041(7) 


0.9314(9) 


1.430(2) 


J1909-3744 


0.898(3) 


0.351(5) 


0.959(4) 


1.464(8) 


B1937+21 


0.862(4) 


0.034(3) 


0.862(4) 


1.451(9) 



Note. — Numbers in parenthesis represent the statistical uncertainty (one standard deviation) in the last digit quoted. All errors are 
normalized by that predicted for conventional pulsar timing based on the total intensity profile alone. 



Table 2 

Conventional Systematic Timing Errors 



Pulsar 


Gr (ns) 


13s 


Is (%) 


5s n 


ref 


J0437-4715 


130 


0.0015 


0.3 


0.2 


5 


J1022+1001 


aeo'^ 


0.011 


2.2 


1.3 


1 


J1713+0747 


180 


0.13 


26 


15 


4 


B1855+09 


530 


0.13 


26 


15 


3 


J1909-3744 


74 


0.015 


3.0 


1.7 


2 


B1937+21 


140 


0.018 


3.6 


2.0 


3 



"Extrapolated to 60min integration. 

References. — (1) Hotan, Bailes, & Ord 2004; (2) Jacoby et al. 2005; (3) Lommen & Backer 2001; (4) Splaver et al. 2005; (5) van Straten 
et al. 2001. 
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W. van Straten 



6k, and ellipticities, e^, of the feed receptors, the differen- 
tial phase, (/), the differential gain ratio, 7, and the absolute 
gain, G, specified in units of an intermediate reference volt- 
age that can later be calibrated to produce absolute flux 
estimates. (Note that in van Straten (2004), 7 refers to the 
boost impact parameter, not the differential gain ratio.) 

In each 500 kHz channel, it is possible to estimate the 
ellipticities and orientations of the feed receptors with an 
uncertainty of only 1°. The apparent correlation between 
those angles indicates that the non-orthogonality of the re- 
ceptors, 5$ and 5^, is better constrained than their absolute 
positions. As a function of frequency, the differential gain 
ratio, 7, exhibits a peak-to-peak variation of about 10%. 
If not properly calibrated, these deviations will result in 
severe systematic distortion of arrival times derived from 
the total intensity profile (cf. Table 2). Furthermore, the 
effect on the frequency-integrated profile may vary ran- 
domly in time due to interstellar scintillation of the pulsar 
signal. Figure 2 emphasizes the importance of calibrating 
with sufficiently high frequency resolution to avoid band- 
width depolarization (van Straten 2002), an irreversible 
effect that cannot be modeled in the current formulation 
of matrix template matching. 

6. CONCLUSIONS 

The future of high-precision pulsar timing is inextrica- 
bly linked with advances in high-fidelity polarimetry. For 
one out of the six selected pulsars, the improvement in 
precision yielded by matrix template matching is better 
than that produced by doubling the integration length or 
instrumental bandwidth. In the analysis of the conven- 
tional method of scalar template matching, conservative 
levels of instrumental distortion are predicted to produce 
systematic timing errors of the same order as the rms tim- 
ing residuals in the current best data sets. These errors 
are completely eliminated by the matrix template match- 
ing method. Therefore, it is expected this technique will 
perform better than current conventional methods in the 
majority of experiments. Furthermore, the method pro- 
vides a new means of fully calibrating the instrumental 



response using a single observation of a well-determined 
pulsar. 

Implicit in the application of the matrix template 
matching method is the assumption that the average po- 
larization intrinsic to the pulsar does not vary significantly 
over any timescale of interest. Although a variety of stud- 
ies have investigated the large fluctuations in the polariza- 
tion of single pulses (e.g. Edwards (2004), Karastergiou & 
Johnston (2004), and references therein) no research on the 
long term properties of pulsar polarization has been pub- 
lished to date. Most likely, any process that effects phase- 
dependent changes in average polarization would similarly 
alter the total intensity profile. Therefore, a systematic 
study of the long term stability of millisecond pulsar po- 
larization would make a valuable contribution toward the 
major science goals of this field, such as the detection of 
low-frequency gravitational radiation and the verification 
of relativistic gravity in the strong-field limit (Cordes et al. 
2004; Hobbs 2005, and references therein). 

The impact of polarization on pulsar timing will become 
increasingly apparent in a larger number of experiments as 
instruments with greater sensitivity are employed, and it 
is imperative to incorporate a sophisticated treatment of 
polarization in high-precision timing analyses. All of the 
software required to perform matrix template matching is 
freely available as part of PSRCHIVE (Hotan, van Straten, 
& Manchester 2004), which has been openly developed in 
an effort to facdlitate the exchange of pulsar astronomical 
data between observatories and research groups. 

I thank F. Jenet for contributing valuable advice and 
feedback throughout this research, and gratefully acknowl- 
edge the thoughtful criticism of the referee. Polarization 
data were provided by S. Ord, I. Stairs, and the data base 
of published pulse profiles maintained by the European 
Pulsar Network. The Parkes Observatory is part of the 
Australia Telescope which is funded by the Commonwealth 
of Australia for operation as a National Facility managed 
by CSIRO. This work was supported by NASA under grant 
NAG5-13396. 
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Fig. 2. — Best-fit instrumental parameters as a function of observing frequency. From top to bottom are plotted the orientations, and 
cUipticities, e^, of the feed receptors, the differential phase, <p, the differential gain ratio, 7, and the absolute gain, G. In the top two panels, 
black and gray correspond to receptors ro and ri, respectively. 



